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Abstract 

The variational properties of the scalar so-called "Universal" equations are reviewed 
and generalised. In particular, we note that contrary to earlier claims, each member of 
the Euler hierarchy may have an explicit field dependence. The Euler hierarchy itself is 
given a new interpretation in terms of the formal complex of variational calculus, and 
is shown to be related to the algebra of distinguished symmetries of the first source 
form. 







1 Introduction 



The Universal Field models proposed in ^, ^, ^ are a class of integrable field theories 
with a wide variety of attractive features: 

• they may be formulated in an arbitrary number of space-time dimensions; 

• they are either diffeomorphism or reparametrisation invariant in the dependent vari- 
ables; 

• they are Lorentz and Euclidean invariant (and in fact G'L(m,R) invariant) in the in- 
dependent variables; 

• they are derived from an infinite number of inequivalent Lagrangians. 

Furthermore, the scalar versions of the theory have been shown to be linearisable. One of 
the scalar theories is a direct generalisation of the Bateman equation [^] and is linearisable 
by a Legendre transformation in the manner described in [Q. 

The business at hand is to describe and explain the strange variational properties of 
the scalar Universal equations. Consider a scalar field (p dependent on m space(-time) 
coordinates. The variational derivation of the Universal equations presented in |Q] is based 
on the idea of the generic Euler hierarchy. Consider a Lagrangian density which only 
depends on the first derivatives of 0. Compute the variation of this and construct a new 
function, 



where 8 is the Euler variation and Fi is some real-valued function dependent only on the first 
derivatives of (j). Then consider £i to be the Lagrangian for a new equation, £Ci. Continue 
the process to arrive at: 



This sequence terminates at the (m + l)th stage: SCm = 0. Furthermore, at the mth stage 
the resulting Euler-Lagrange expression factorises, and one of the factors is "universal" — 
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independent of the details of the initial and intermediate Lagrangians. On setting this Euler- 
Lagrange form to zero, we arrive at an equation of motion equivalent to the Monge- Ampere 
expression, 

AM-A = detH = 0, (1.3) 

where H is the Hessian matrix of the dependent variable 0: Hij = (pij. 

The most interesting special cases of this construction are referred to in [Q] as generalised 
Bateman equations. The idea is to use an initial Lagrangian £o which is homogeneous of 
degree one in the first derivatives of 0, and to restrict all the to have this property too. 
Performing the procedure described above leads to the cessation of the hierarchy a stage 
earher than described above: £Cm~i = 0. Remarkably, the penultimate Euler-Lagrange 
expression EC,m-2 is again a product, and again one of the its factors is "universal". Setting 
this expression to zero gives the Universal Bateman Equation of []T| and sequels. 

The equation can be written. 



A = det 



/ 0^, • • • 0x™ \ 

0X1 0X1X1 ■ ■ ■ (pXiXm. 



0, (1.4) 



or in components, 

A = £n...im^jl...j™0ii0ii0i2i2 • • • = 0- (1-5) 

The component form makes it easy to verify that A can also be expressed as: 

A = tT[GH^), (1.6) 

where the matrix G has components 

Gij = (j)i(pj (1.7) 

and 

H^ = adi{H) (1.8) 
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is the classical adjoint matrix of H. Alternatively, we could define a new matrix U such that 
(|L6|) is equivalent to the equation: 

/\ = ii{UH). (1.9) 
The explicit form of U is easily deduced from (|1.5|) . It is: 

U = £h...i^ej^,„j^(f)i^(f)j^(j)i^j^ . . . (1.10) 

A supplementary assumption is that det H is non-vanishing. 

Our mission is to try to understand and generalise the generic and Bateman hierarchies 
described above from the point of view of the standard theory of variational calculus as 
presented in, for example, the book by Olver 0. In particular, we wish to know under what 
circumstances we can introduce dependence on the field (p into the initial Lagrangian and 
the multiplier functions F^. We will find that only in certain circumstances does the Euler 
hierarchy define a "universal" theory. In addition, we will think about what the iterative 
procedure means in terms of the formal variational complex, described in or the work of 
Anderson [0. We will find that the iterated equations of motion are related to determining 
equations for distinguished symmetries of higher members of the hierarchy, and we will justify 
this connection with a calculation of the generalised symmetries of the Universal equations. 

Basic references (apart from the original papers [T|, ^ ^, ^) are Olver the review 
articles by Anderson [0, |[ and his recent work with Pohjanpelto |T0|- Another view of the 



universal equations and their variational symmetry properties is given by Grigore [O, O . 



2 Lagrangian Properties 

We begin by analysing the Euler hierarchy in the language of variational calculus presented 
in . Initially, we will only assume that the initial Lagrangian depends on the field and its 



first derivatives. In the language of jet bundles |]T3|, we ask that the Lagrangian is a smooth 
function Cq : J^vr -^W, where J^tt is the first jet bundle of the bundle tt. In this case we 
are restricting our attention to the trivial bundle tt :R'"+^ -h>M™. We make no assumptions 
about any symmetries of the resulting action. This language will be useful later. 



From Co, we derive an Euler-Lagrange form: 

A, = Dl^il)=SU (2.1) 

From this we construct the next Lagrangian, 

A = FiAi, (2.2) 

where we assume that Fi depends only on (p and its first derivatives. The next source form 
follows readily using the product formula: 

£(^i^2) = D>^(^2) + D>^(^i), (2.3) 

where D* indicates the adjoint of the Frechet derivative 0. In this case: 

A2 = £:£i = D1^(Fi) + D:^^(Ai). (2.4) 

Now, the Helmholtz condition of the calculus of variations states that an equation is an exact 
Euler-Lagrange form if and only if its Frechet derivative is self-adjoint. Hence: 



DA,(i^l) = DA.(Fi) 

dAi 



The second term on the right-hand side of ( |2.4| ) just turns out to be 

dFi ^ fdFi 



Dj{Fi) (2.5) 



DK(A0-A,^-A^AJ, (2.6) 



since Fi only depends on (p and its first derivatives. 



Putting (p^ ) and (|2.6|) into ( p.4|) , a short calculation tells us that: 



The important point to notice is that this construction guarantees that there are no deriva- 
tives of order higher than two in the resulting equation of motion. This means that when 



the process is continued, the analysis is essentially unchanged. We are led to the iterative 
formula for the {k + l)th Euler-Lagrange form: 

, ^ ^ dAk f ^ dF.\ dAk dFk , 
^ dcj) \ d(pi ^ ) d(j)i d(p ^ 

There is no such recursive definition if Fk depends on second or higher derivatives. 

We are interested in sequences generated by this kind of recursion which terminate after a 
finite number of iterations. The expression ( p.8|) simplifies greatly if we restrict attention to 
Fk that are (1) independent of (/>, and (2) homogeneous of degree one in the first derivatives. 
Then we find: 

A,,..A,£f, + ^B,(g^„ ,2.9) 

or more symmetrically, 

^k+i = ■w—(Pik(Pji - Ak(pij . (2.10) 

\d(f)ki J d<j)id(j)j 

This is precisely the recurrence found by Fairlie and Govaerts (0]) in their treatment 

of the generic hierarchy. They found that if Cq and all the F^ were independent of (p then 

the sequence terminated in a Monge-Ampere equation. This sequence only produces such 

a result if it is assumed that the initial Lagrangian is independent of <p, in which case the 

condition (2) can be relaxed anyway, so there is no new information. Otherwise, ( |2.9| ) 

fails to reproduce the generic hierarchy. Termination of the hierarchy seems to depend on 

certain linear-algebraic properties induced by the fact that the second derivatives appear 

homogeneously at each stage of the sequence . This does not happen if the F^ depend on 

0. So far, all attempts to generalise the construction of the generic hierarchy by relaxing this 

requirement have failed. Sample computer calculations carried out on MAPLE show that 

such constructions do not terminate in universal, exact source forms at either the (m — l)th 

or mth stages. 

We proceed to generate the equations of the generic hierarchy under the restriction that 
the Co and Fk are 0-independent. We know that the first equation ( |2.1| ) has the explicit 



form, 

which can be written: 

Ai = -tr (HMo) , (2.12) 

where Mq is the Hessian matrix of Co with respect to its dependence on first derivatives of 
(j). It is then a straightforward matter to apply the recursion relation ( p.lO|) to this starting 
term. Using the notation. 



^3 

Pk = HMk, (2.13) 



the first few terms of the hierarchy are: 
Ai = -tr(Po), 

A2 = tr(Po)tr(Pi)-tr(PoA), 



A3 = -tr (P2 (PoPi + PiPo)) + tr (Po) tr (P1P2) + tr (Pi) tr (P0P2) 

+ tr (P2) tr (PoPi) - tr (Pq) tr (Pi) tr (P2) , 
A4 = -tr (P3 (P1P2P0 + P0P1P2 + P2P0P1 + ^0^2^! + P1P0P2 + P2P1P0)) 
+ tr (P1P2) tr (P0P3) + tr (P2P3) tr (PqPi) + tr (P1P3) tr (P0P2) 
+ tr (Po) (tr (P1P2P3) + tr (P3P2P1)) + tr (Pi) (tr (P0P2P3) + tr (P3P2P0)) 
+ tr (P2) (tr (P0P1P3) + tr (P3P1P0)) 

- tr (Po) tr (Pi) tr (P2P3) - tr (Pq) tr (P2) tr (P1P3) - tr (Pi) tr (P2) tr (P0P3) 
-tr(P3)A3 

; (2.14) 

The subsequent members of the hierarchy become combinatorially more complicated. 

We can use the recursion relation as the basis for a plausibility argument that demon- 
strates the vanishing of the A^ at a certain stage, and hence when the Fk is the characteristic 
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of a conservation law. We wish to know when vanishes for any possible choices of Fk- 
(Of course, we are only interested in those cases when the Hessian of Fk is non- vanishing.) 
This amounts to solving the matrix differential equation, 

HT^'H = AkH, (2.15) 

where, 

The equation (|2.15| ) is solved by any of the form A = (some factor) x det H, where "some 
factor" is independent of the second derivatives of (p. We know from the variational calcu- 
lation that the second derivatives only enter the source forms such that Afc is homogeneous 
of degree k in the second derivatives. Therefore, this solution can only work at the mth 
stage. We will see in the next section that any function of the first derivatives of (j) is the 
characteristic of a generalised symmetry of Am-a- It follows immediately that taking an 
Euler variation of Aj\/_^ as described above will give a zero result. 

A more rigorous proof of these results is presented in using an explicit calculation. 

This completes the discussion of the generic case. Now we turn to the Bateman case. 
The reader should keep in mind Euler's theorem on homogeneous functions, which states 
that, if ^(x) is homogeneous of degree a in its arguments x = (xi, ...,x„), then: 

In our case, we assume that £o is homogeneous of degree one in the first derivatives of 0, 
meaning that, 

<l>^Jr^ = U (2.18) 

It is easy to derive further homogeneity properties of the Euler variation and the derivatives 
of Cq. 

In particular, the first derivative of Cq with respect to (pi is homogeneous of degree zero: 



'd(j)id(pj 



0, (2.19) 



which imphes the singularity condition: 



det(Mo) = 0, (2.20) 

since <p is arbitrary. 

So we proceed to discussion of the Bateman hierarchy. Now we will assume that the initial 
Lagrangian and the Fk have arbitrary dependence on (p and are homogeneous of degree one in 
first derivatives of (f). This greatly civilises the shapes of both Ai and A2 from their original 
forms ( |2.1| ) and ( p.7[ ). If we use the various properties stemming from Euler's theorem, we 
recalculate Ai and A2 to be, 



((9 A \ f)'^ F 

d<f)M J d(f)id(f)j 



from which we deduce precisely the same recurrence ( |2.1Q|) without the restriction that the 
Fk need be independent of (j). 



The recursive procedure now defines a set of equations identical in form to (p.l4|) , and so 



the proof of Fairlie and Govaerts [Q is still valid. The equation governing the termination 
of the sequence ( p.l5| ) is satisfied by (|1.5|). 



There is a nice interpretation of all this in terms of the theory of the Euler-Lagrange 
complex studied in [0. Let Q* denote the algebra of differential forms on the infinite jet 
bundle J°°7r. This algebra is bigraded into, 

= 0^]('-'^), (2.22) 

r,s 

where the spaces il*^^'*) contain forms with the local coordinate expression: 

f[(p] dx'' A ... A dx'"- AOf^ A...A 

the forms 6 being the various contact one-forms. The various differentials and subcomplexes 
of this space are discussed in [0. Alternatively, see the book by Saunders |jl3|. 
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The starting Lagrangian is a form Aq G ^1^"^'^^ The corresponding variation is the source 
form Aidcf) A u = S{X) G Q"^'^. Now consider taking a Lie derivative (denoted by L) of Ai 
with respect to the prolongation of some vector field Xi whose evolutionary form is Xp-^, 
with characteristic Fi. For an arbitrary equation A, we have the Cartan formula: 

LprXi(A d(f)ALo) = Sv {Xp^j A d(t)Auj) + Xp^j 6v (A d(/)Auj). (2.23) 

The differentials 6v are those of the Euler-Lagrange complex defined in or 0. A proof 
of ( p.23| ) can be found in |p. Each term has a simple interpretation in the calculus of 



variations. The left hand side vanishes identically if Xi is a distinguished symmetry of A. 
The first term on the right is an Euler variation, and so it vanishes if Fi is the characteristic 
of a conservation law. Finally, the 6v in the third term is just a Helmholtz operator, so it 
vanishes if A is an exact Euler-Lagrange form. Not surprisingly, this is the central formula in 
the recent studies of the generalised Noether theorem by Anderson and Pohjanpelto P, p!0|. 



Returning to our example, ( |2.23| ) takes on a particularly simple form when applied to 
Ai. Since Ai is an exact Euler-Lagrange form, the Helmholtz term vanishes and we are left 
with the form A2 as defined in (|2.4| ). This interpretation holds for all A^, so we can view 
the Euler hierarchy as repeated application of the Lie derivative to successive source forms. 
So the recursive definition becomes: 

Afe+i d(t)Auj = Lprv^^ (Afc d(t)Au). (2.24) 

The "universal" theory rests on the observation that at a certain stage the all source forms 
so defined, in either the generic or Bateman hierarchy, are equivalent to one another and 
their flows defined by the Lie derivative are identically zero. In such a situation, the iteration 
vanishes identically, and this yields a product expansion of the type (^) which is equal to 



zero. This indicates that Ck is a total divergence, and hence F^ must be the characteristic of 
a conservation law. Explicitly, we deduce that any function of the first derivatives of is a 
generalised symmetry of the Monge- Ampere equation and that any function of homogeneous 
of weight one in the 0j but with arbitrary dependence on is a generalised symmetry of the 
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Bateman-type Universal equation ( |1.4| ). These results will be proved explicitly in the next 
section. 

Changing emphasis, we can view the Fk as unknown functions to be determined, and 
then we can interpret A2 as the equation determining the distinguished symmetry algebra 
of the Euler-Lagrange form Ai. The third and subsequent source forms A3, ... , A^ are a 
set of nested equations determining distinguished symmetry algebras for their immediate 
predecessors. It would be interesting to know what information, if any, can be gleaned from 
these higher equations about the symmetries of the original equation Ai. 

3 Generalised Symmetries 

In order to confirm our understanding of how the generic and Bateman hierarchies terminate, 
it will be helpful to know more about the generalised symmetries of the equations (|1.3D and 



( |1.4|) . To that end we will carry out a detailed analysis of each of each equation in the 
manner described in Olver's book p. We will look for first-order generalised symmetries 
— in other words symmetries whose evolutionary characteristic depends on the Xi, cj) and 
the first derivatives of (f). The reader is reminded of the standard expansions of the total 
derivatives D., and Dj 



^ ^ OF dF ^ 
oxi dcpj 

d'ip d^F OF d^F d^F 

^ii^ = ^ a a - (t>Kj + -W—<t>Jij + o o , 0J» + o , o , 0J»0i^j^ (3-1) 

OXiOXj dXi04>K Ocpj OXjOCpj OCpjOCpK 

where the multi-indices J, K will actually only have length zero or one in the cases we are 
considering. 

We will begin with the equation It turns out that the first-order generalised 

symmetries of this equation span a rather large infinity of possibilities. Using the standard 
procedure, we seek a generalised symmetry in the evolutionary form, 

VQ = q|-. (3.2) 
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The symmetry equation is quite straightforward: 

prvQ(AM-A) = ^^A,QM = 4A,-Q[0]. (3.3) 

The symbol pr denotes the infinite prolongation. On expanding the total derivative, we find 
that the coefficients of the third derivatives of <!) vanish on solutions of ( |1.3|) (after taking into 
account the first prolongation of the equation of motion) and the first symmetry constraint 
comes from the term of order (m + 1) in the second derivatives, which gives: 

Hlj^§-<l>,kc^,i = 0, (3.4) 
d(f)kd(pi 



on solutions. Writing, 



we find that ( ^.41 ) becomes: 

tr (h^HHqH) = det H tr {HqH) = 0, (3.6) 

and this is identically true on solutions of Am- a = 0. This pattern continues for the terms 
of order m in second derivatives, and we find that a Q with arbitrary dependence on the 
first derivatives is a generalised symmetry for (|1.3| ). This is the expected result given the 
interpretation of the generic hierarchy and confirms its termination at the (m + l)th stage. 

Now we turn our attention to the symmetry algebra of the Universal equation (|1.4|) , 
which turns out to be a little more restrictive than in the generic case. The determining 
equation is, 

|^(AQ[0]) + ^(A,gM) = 0, (3.7) 

on solutions of A = 0. Again, we will only look for first-order generalised symmetries. 
This should be sufficient to confirm our interpretation of the Bateman hierarchy in the 
previous section. Furthermore, experience with the two dimensional case suggests that 
the symmetries dependent on first derivatives are responsible for the linearisability property. 
Given this assumption, (|3.7|) expands to become. 



^h-i^^h-j,AiAh<Pi2j2 ■ ■ ■ <Pimjm + ("1 - l)U kl5 (^ij)(kl)D ijQ = 0, (3. 
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where, 

Ukl = £h...im-lk£jl...jm.ll(Pil4'jl(pi2j2 ■ ■ ■ 0i,„-lim-l (3-9) 

is the sort of matrix appearing in 

Following the usual algorithm, we must set the left side of (|3.8| ) to zero order-by-order 
in the derivatives, taking into account the equation A = 0. Using the first prolongation of 
A, we find that third order derivatives automatically vanish, and so the first task is to find 
the condition for the vanishing of terms of order m in the second derivatives of (p. 

Extracting the relevant term from ( p.8|) , we define a matrix S such that: 



0, (3.10) 



on solutions of A = 0. To simplify this, we use the form ( |1.9| ) to rewrite A using the cyclicity 
and linearity of trace and the properties of the classical adjoint: 



m det H ^ ' 



A = — ^ ^T^trfSi/^). (3.12) 



and then, by associativity. 



m(m — 1) det H 

On comparison with ( |1.6| ) we find that S = m{m — 1) det H G and the "on-shell" symmetry 
condition is satisfied by any S =(some scalar factor) x det H G. This produces a symmetry 
condition, 

0.'^,7rr^ = 0, (3.13) 
d4)id(t)j 

which is satisfied by any Q homogeneous of degree zero or one in the first derivatives of 0. 
This confirms the facts deduced earlier from the Bateman hierarchy. 

According to a theorem of Kumei and Bluman |]l4l, the existence of a generalised symme- 



try which depends on the solution to a linear equation such as (|3.13|) guarantees the existence 
of a linearising transform for (|1.4| ) . This is indeed the demonstrated in [ffl. 
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Having disposed of the term with m second derivatives, we need to equate the term with 
m — 1 second derivatives to zero. The relevant equation can be written as: 

4||</>, + 2(m - [-i;§fj^n + ^J^<P^r}j = 0. (3.14) 

This can be solved by a characteristic of the form, 

Q = g{<P)F{r]), (3.15) 

where g is an arbitrary smooth, real- valued function, rj = XiCpi (no sum) and F is a smooth, 
real-valued function which respects the homogeneity properties that we have decided on for 
Q. Notice that the diffeomorphism symmetry of the equation is included in this solution. 

Finally, to get rid of the remaining terms, we use ( 3.15|) and arrive eventually at the 
condition: 

£ii...i™ei,...i^(/>n(/>ji0i2i2 ■ ■ ■ Mj ■ ■ ■ (t>imjm F" + 2g' F' + g"F) = 0. (3.16) 
This is identically true due to the antisymmetry of the e symbol. 



4 Conclusion 

The generic Euler hierarchy and the special case of the Bateman hierarchy have been inter- 
preted as a sequence of iterated Lie derivatives for the distinguished symmetries of a large 
class of Lagrangians. It has been shown that the Bateman hierarchy alone admits explicit 
dependency on the field due to its homogeneity properties, but the form of the generic hierar- 
chy is unchanged from the analysis in . The termination of both hierarchies is guaranteed 
by the generalised symmetries of their associated universal equations, which in turn imply 
the existence of an infinite class of conservation laws for each equation. 

Many questions remain unanswered about the geometrical meaning of the iteration of 
the Lie derivatives, and the presence of any helpful algebraic properties. It also remains 
to apply the procedure to more general geometrical constructions than the trivial bundle 
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TT considered here. It is noted that other equations follow from similar procedures, namely 
certain multidimensional generalisations of the Born-Infeld equation ||15| . 
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